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ABSTRAC 



An experiment using the focusing properties of a nonuniform 
axially symmetric magnetic field to study large angle (0 > 32°) 
scattering of ions from atoms and molecules is described. The 
scattering cell is placed on the magnetic field axis and all ions 
in the momentum interval j Apl at |p| scattered into a conical shell 
AO at 0 are counted by a detector placed further down the axis. The 
solid angle is increased by a factor of several hundred over con- 
ventional scattering techniques of comparable angular resolution. 
The magnetic vector potential and the trajectory equation for a 
charged particle in the field are derived. The computer program 
which solves the trajectory equation is given along with a detailed 
description of the apparatus. Our measurements of the absolute 
differential scattering cross sections for Li^ on He from 36° to 
41 *7 at 200 eV and 3 OO eV and for He on He from 40° to 54° at 
300 eV and 400 eV demonstrate the feasibility of the technique. 
Using the measured cross sections for Li^ on He and the computa- 
tional method of Firsov we have calculated the interatomic poten- 
tial function for nuclear separations of O. 15 X to l.lS. 
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CHAPTER I 



INTRODUCTION 

In recent years the techniques for measuring differential 
scattering cross sections for various atomic and molecular pro- 
cesses have been improved to the point where very precise mea- 
surements can be made in an almost routine manner. In the con- 
ventional approach a target of N^ particles is confined within a 
surrounding vacuum. A well collimated beam of nearly monoener - 
getic ions of current I^ is projected into the target. A de- 
tector, subtending a solid angle dQ (0) = Sin 0 dO 60 at the 
interaction region, is placed at an angle 0 to the incident ion 
beam and the number of particles entering the detector,- N^, is 
measured. The scattered particles are usually detected as a 
current where I^ = N^dn(O)^ From these measured quantities, the 
differential scattering cross section in the laboratory coordinate 
system can be calculated as 

da(0) = I /l.N^dQ(O). 

' ' sit 

However, in order to obtain reasonably good angular resolution, 

the solid angle includes only about 1^ of the available 360 ° of 

azimuthal angle 0 into which ions scatter. This restriction re- 

-Zf 

suits in fairly small solid angles of about 10 steradians, 

-20 2 

which sets a lower limit of around 10 cm on observable labora- 
tory differential scattering cross sections. By modifying the 
techniques of thin magnetic lens p-ray spectroscopy, we have con- 
structed an experimental apparatus which collects essentially all 
ions scattered into a 360 cone dO wide at 0. The solid angle in 
this apparatus is about 300 times larger than that of other methods. 
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allowing us to observe differential scattering at laboratory 
angles greater than 32 degrees and cross sections substantially 
smaller than any previously measured « 

The apparatus p shown schematically in Figure 1, was devel- 
oped to make use of the focusing properties of an axially sym- 
metric nonuniform magnetic field. The scattering cell exit slit 
o 

is open 290 of the aximuthal angle, allowing ions of a given p 
and 0 but scattered at most any azimuthal angle to be focused 
at the detector by the magnetic field. Both the magnet current 
and the distance between the scattering cell and detector are 
Variable. Various combinations of the magnet current and the 
detector distance make a wide range of reactions observable with 
the apparatus. The details of the apparatus are discussed in 
Chapter III. 

The differential equation for the trajectory of a charged 
particle in a nonuniform axially symmetric magnetic field and 
the magnetic vector potential for this field are developed in 
Chapter II. The trajectory equation is numerically integrated 
to give the values of the magnetic field and detector position 
necessary to observe charged particles scattered at various angles 
and energies. The solid angle is calculated from the computed 
trajectories and is a function of the angle of scatter, the en- 
ergy of the scattered ion, the magnetic field, and the detector 
distance. This calculation is discussed in detail in Chapter II. 

• • t' 

The elastic differential scattering cross sections for He 
on He at 300 eV from ^0° to 52°, and for He^ on He at 400 eV from 
42° to 53° were measured using this method and the results are in 
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ION SOURCE 
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Figure 1. Schematic of Experimental Apparatus 



essential agreement with measurements of Lorents and Aberth^ who 
used a standard technique. The elastic differential scattering 

6 4 * 

cross sections measured for Lx on He at 200 eV and 300 eV from 
34^ to 4l.7^ also compared favorably to extrapolated values of 

2 

cross sections measured at smaller angles by other investigators . 

These experimental results are given in Chapter IV. 

Finally, the elastic differential scattering cross section 

6 + 

for Li on He measured by our method, along wxth data taken at 
smaller angles by other workers, was used to calculate the ion- 

atom potential function for this system using a method developed 

3 . . 

by Firsov . This calculation is discussed in Chapter V. 

Based on these results we conclude that the feasibility of 
the technique has been established and its usefulness in obtaining 
new data has been demonstrated. 
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CHAPTER II 



KINEMATICS OF CHARGED PARTICLES IN A NONUNIFORM 
AXIALLY SYMMETRIC MAGNETIC FIELD 
In this chapter the equation of motion of a charged particle 
in a uniform magnetic field is solved to show the explicit re- 
lation among the following three variables: the vector momentum 
—* — ♦ 

of the ion, p, the magnetic field, B, and z^, the distance from 
the origin at which a detector must be placed to observe the ion. 
The expression for the solid angle is obtained in closed form 
and its dependence on the various parameters is shown. 

The orbit equation for the axially symmetric nonu niform 
magnetic field is then derived along with the expression for the 
magnetic vector potential, which is part of the orbit equation. 
Expressions for the angular and energy resolution are found and 
these results are used to find the solid angle. The dependence 
of the solid angle on the various parameters is shown and dis- 
cussed. 

Finally, an experimental method of simulating a particle 
trajectory in a magnetic field using a flexible, current carrying 
wire is described. The results of these measurements were used 
for preliminary design considerations and are compared to the 
orbits which were obtained by numerical integration of the tra- 
jectory equation. 

1. MOTION OF A CHARGED PARTICLE IN A UNIFORM MAGNETIC FIELD 

A longitudinal type homogeneous magnetic field spectrometer 

has commonly been used in the analysis of p spectra in nuclear 
A 5 

physics ^ . The theory of motion of a charged particle in a 
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uniform magnetic field was originally published by Busch in 1926, 

upon which the technique was based. 

The trajectory of a charged particle in a uniform field is a 

helix, where the particle, after one complete revolution of period 

T, will focus again on the axis from which it originated. As seen 

in Figure 2, the distance at which the ion returns to the Z 

axis, depends on the momentum, p, the angle, Q, and the magnetic 
— ♦ 

field, B. The relationship between these quantities can be derived 
as follows: 

Let p = vector momentum of an ion at the origin 

of the coordinate system at t = 0 

and q = charge of ion c 

The magnetic field is described by 

A A A 

B = B r+B, ^+B Z 
r 0 z 

where B = B, = 0, B = constant, 

r ‘ 0 ^ z 

The boundary conditions seen from Figure 2 are; 



at 


o 

11 


K 

II 

o 


z = 0, p^ = p cos 0, p^ = p sin 0, 


and at 


0 0 

H 

II 


r = 0, 


z = Z . 



In order to get the relation Z^ = f(p,0,B) we must integrate 

— » — * 

the force equation, F = q v x B, and apply the boundary conditions 
at t = 0, and t = T. Thus 
F = q V X B, 
and in component form 



F 

r 



= m 





^ 0 Z 




( 1 ) 



l6 



DETECTOR 
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Figure 2. Orbit of a Charged Particle in a Uniform Magnetic Field 



F - = m 
0 



.(i)(s> .01 



q V B = 
^ r z 






d^z „ 

F = m ;r = 0 . 

^ dt^ 



( 3 ) 



Integrating ( 3 ) and applying the boundary conditions at t 
gives 

m -tt - constant = p cos 0 
dt 



= 0 



j^(p/m) cos oj t. 



( 4 ) 



From (2) , upon multiplying by r and integrating 



m 



2r 









m 



dt 



dt 



d , 2, 
— dt <' ) 



dt 2m 



= constant = k. 



( 5 ) 



Now using ( 5 ) in (1) gives 



d^r 2 

m - mrk = ~qrk 



= 2rmk 



dt 



d^r ■ 2 

— - + rk =0 

dt 



which has the solution 



^ ikt ^ ■ -ikt 

r = e ^2 ^ 



From the initial conditions at t = 0 



= -C^ = p s in 0/2ikm 
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giving 



r 



p sin Q 
km 



s in kt . 



( 6 ) 



Now eliminating t from equations ( 4 ) and ( 6 ) gives 



r 




and at t = T, this reduces to 




_ p cos 0 
2n ^ 

q 



( 7 ) 



Therefore, a detector at Z with a field B will ”see" a particle 

o z ^ 

of charge q scattered with momentum p at an angle 0 with the z 
axis. Notice that there are two ways to observe different 



detector, (2) fix the detector at a particular value of Z and 



a detector to collect all ions within a given momentum interval 



make observations, then, one would simply choose the values of 

p cos 0 to be observed, solve the equation for various values of 

Z and B that will allow them to be seen, and set the detector 
o z 

and field accordingly. 

The solid angle in a spectrometer having asrimuthal symmetry 
depends on 0 and p as follows: recall that dQ = 2n sin 0 dO 

where we have integrated over d0. An expression for dO can be 
obtained by differentiating equation (7) as follows: 



p cos 0 values: (1) hold B^ constant and change Z^ by moving the 



o 



vary B^. Thus a knowledge of equation ( 7 ) allows one to position 



I Ap| at 1 p 1 scattered into a conical shell AO at 0. In order to 




and since B^ = constant, then 
2 ^ 




qB, 
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and 



dO = 



dp 



qB dZ 
z o 



p tan 0 2tt p sin Q 



However 



m 



p = (2mE) or dp = (— ) dE 



giving 



dO = 



dE 



qB dZ 
z o 
U 



2E tan 0 2n(2mE)'^sin 0 



Using this the solid angle is expressed as 

/ 2 \^ 

, ^ cos 0 V ^ z \ 

< — i — > 



( 8 ) 



Therefore, as seen from equation (8) the solid angle depends on 
the angle of scatter ^ 0, the trajectory intercept spread, dz^, 
allowed by the detector, and the energy of the scattered ion, E. 

As an example, consider the elastic collision He"*" on He 
where 0 = 30°, E = 300 eV, dE = 5 eV, = 200 gauss, and 
dZ^ = - 3inm* Then dQ = O.O 58 steradians, a value typical of ' 
solid angles for the nonu niform magnetic field apparatus actually 
used . 



2. ANALOGUE MEASUREMENT OF ORBITS 
Prior to obtaining any numerical solutions for the orbit 
equation, the technique of simulating the orbit of a charged 
particle in a magnetic field using a current carrying flexible 
wire 'was employed to obtain preliminary design information on the 
various experimental parameters • 

7 

The experiment, carried out by Kelly , is based on the fact 
that a current carrying wire consists of charged particles moving 
with some velocity along the wire, causing the wire to experience 

a V X B force in a magnetic field, the same type of force charged 
par t ides encounter . 
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Consider a small length, ds , of flexible wire carrying a 



current, i, and under a tension, T, in a magnetic field, B. 




p = particle momentum 
q = particle charge 
p = radius of curvature 
m = mass/length 
i = current 
T = tension 

— I 

T = unit vector along 
tangent 

T) = unit vector along normal 



Then since the system is in equilibrium the sum of the forces is 
zero, giving 

— » — > — » — ♦ 
ids(T X B) - T T| - mg ds = 0 

from which the curvature of the wire is shown to be 



in 

p 




X B) 



mg 

T 



( 1 ) 



To get the path of a charged particle in a magnetic field, recall 
that 

^ = 



dt 



q(v X B) 



p V - 

but T = ^ — and ds = vdt, for constant velocity, v* 

p V 

Therefore v ^ = vp ^ = qv(r x B) , 



however 



dT ^ T) 
ds P 



and the curvature of path of a charged particle is 



1 = S f T 



P P 



(T X B) 



( 2 ) 
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Therefore the wire and particle paths will coincide when 




( 3 ) 



Now for our case the magnetic force is usually 100 times the 
gravitational force, hence the gravitational term can be neg- 



Thus when equation (4) is satisfied it allows the simulation of 
charged particle trajectories by a current carrying wire under 
tens ion . 

The experimental arrangement used to make the wire orbit 

measurements is shown in Figure 3. An uninsulated wire con- 
sisting of l4 strands of number 44 copper wire weighing 0.24 

grams /meter was used. Currents of up to 4«5 amps could be 

passed through the wire, however at the higher currents the 
wire oxidized and became less flexible. The tension in the 
wire was normally 1000 dynes and was measured by a modified 
ammeter movement which was calibrated against a set of weights. 
This method was chosen over the usual pulley arrangement since 
the frictional effects in the pulley could not be sufficiently 
reduced. The departure angle, 0, was measured to an accuracy of 
^ 30 minutes of arc by reading the longitudal position of a copper 
sleeve at the instant it made electrical contact with the wire. 

The measured departure angle was then corrected to account for the 
curvature of the wire between the teflon ball and the copper 
sleeve . 



lected and equation (3) reduces to 



I = E 
i q 



(4) 
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Figure 3- Schematic of Wire Orbit Method 



For a given magnetic field and distance 2^ the current 
through the wire was gradually increased until the desired de- 
parture angle was obtained. Then the tension and the current were 
measured giving the momentum (for an assumed charge) which the wire 
trajectory represented. The shape of each trajectory was obtained 
by measuring r as a function of Z. Thirty seven different tra- 
jectories were measured and tabulated in reference 7 which also 
includes a detailed discussion of the method. 

A comparison of these measurements with the computed tra- 
jectories showed that the measured maximum r values were about 
17 % larger than the calculated maximum r values and the measured 
Z values were about 8% larger than the calculated Z values* The 
natural stiffness of the wire is believed to be the primary cause 

of the discrepancy. A comparison of the wire orbit and the calcu- 

+ 

lated trajectories is shown in Figure 4 for an ion with a de- 
parture angle of 48.17 degrees and scattered energy of 243 eV. 



3. MAGNETIC VECTOR POTENTIAL 

The trajectory equation has as one of its parameters the 

magnetic vector potential. Therefore, the measured field of the 

magnet must be related to this quantity before the orbits can be 

— ♦ 

found. The well known expansion of A in terms of the field is 

. 6 

found by following the development used originally by Busch . 
Since there are no currents in the region of interest 

V X "b = 0 
and, consequently 

V X (V X A, = 0 
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Figure 4 . Comparison of Wire Orbit and Computer Trajectories for 



is the equation to solve. Now 



V X (V xA) = V(V • A) - V A 

* A |J» ( f I ^ 

and recalling that J - p only, then A = I — ^ 



(1) 



requires that % = A^0 and A^ = A^ = 0. Also, since cylindrical 



symmetry exists, A^ is not a function of 0. Thus in cylindrical 
coordinates, (r,0,z), equation (-1) becomes 



1% + i ^ *0 - „ 

3.= ■ ■ 



Now assume a solution 



then 



A(r,z) ^-0 ^ 



SA ^ 
r — = na r 

Sr n 



-s 2 

S A ^ . n-2 

= Sn(n-l)a r 

^2 ^ ' n 

Sr 



a A ^ ’’ n 

= S a r . 

^2 n 



( 2 ) 



Substituting into equation (2) gives 

n-9 — n-9 '' VI _ 

h S a r =0 



n 



n 




n 


na r 


a r 


n 




n 


n-2 


11 


n 


K 

+ 


u a 


r = 


1 


n 





/ -IX/ -.X ^ ” n 

E (n*-l)(n+l) ar +Sa r =0* 

^ ^ ^ ^ n n 

Each coefficient must be zero and shifting indices in the second 

term gives oo oo 

S „(n+l)(n-l)a r + 

n=0 n n=2 n-2 
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and for 


n = 


0: 


1 

M 

o 

II 

o 

o 

II 

o 




n = 


1: 


O(a^) ~ ~ arbitrary function of z 








!f 




n ^ 


2; 


(n+l)(n-l)a + a . = 0, 
' ' n n-2 



The recurrence relation for the coefficients is 

If 



(n+l)(n-l) 



and since a^ = 0, there are only odd n terms in the solution< 



The general term is 



(-1) 



(2k) 



2k+l 2^'^k*(k+l)! ^ 



This gives 



A(r ,2) = S 



k=0 



^ / -1 N ^ 2k+l . , 

2(-l) ^(2k) 

kl(k+l)* ^2^ 1 



or 



3 5 

" ( TV \ r 

A(r,2) = a^r - ^ + ... 



Now to evaluate the coefficient aj^(2) recall from B = ^ x A 



and the field symmetry that 

1 B(rA) _ 1 a r. -2 " r^ , (IV) r^ 

B -a, -5- + a; iTtTT- ***] 

zrBr r 9 rl lo 1 192 

= P, ^ + ... 

^1 “ ^1 ^ ° 1 192 ’ 

but B (r = 0 ) = 2 a^ giving a (z) = -i- B (r = 0 ). 
z 1 j- d z 



This gives 



A^(r,z) 



= L 



(- 1 ) 



2n+l 3 ^"b (r=0) 

(I) 



n=0 n!(n+l) I ^2' 



2n 



( 3 ) 
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Therefoxe, by knowing B^(Z)^_q we also know ,z) , the quantity 

used to solve the orbit equation. The actual measurements were 
fitted to a 12^^ order polynomial^^ in z such that 

13 , 

B^(z) = S (4) 

With this expression for B^(z )5 the field values were repro- 
ducible to within 0.1% of the experimentally measured values. 

Now using these expressions for and B^(z) in the 

computer program one can calculate the trajectory of a particular 
p/q value for a selected magnetic field. 

4. ORBIT EQUATION IN A NONUNIFORM FIELD 
The thin lens spectrometer with its axially symmetric non - 
uniform field was chosen for this apparatus. It more easily 
produces the large field required to bend ions and it gives more 
spread in the intercept distances, allowing greater resolution. 

The equation of motion for a charged particle in an axially sym- 
metric nonu niform magnetic field and the subsequent employment 

of this field geometry as a p-ray spectrometer has been discussed 

8 9 

by several authors ^ . Our derivation of the trajectory equation 

g 

essentially follows that of Busch , and was first applied to our 
particular problem by Gagliano^^. 

An axially symmetric magnetic field can be represented by 
the curl of the magnetic vector potential having only an azimuthal 
component, A^ (see previous section). The motion of a charged 
particle in such a field is governed by the equation 

— V — ♦ — » 

F=qvxB=qv x(V x A). 
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In cylindrical coordinates, (r, 0, z), the force components are 



rd^r dp 

Fr = n[— - r(^) 1 = q 
dt 



(1) 



*^0 = ™ 



f' 



. dt 



+r 



d^0 

dt^ 



= -q 



fdz\r^0 ' 

, dt / a 2 



1 dr 



+ — 
r dt 



ar 



( 2 ) 



d^z d0 ^ “^0 

= m ^ = q r — 

dt 



(3) 



Now multiplying equation (2) by r and using the relation 



d 

dt 



dz h 
dt Sz 



dr 5 
dt 9r 



(4) 



gives 



d 

dt 



r™r2 d0-, 



q [ 



dt 



1 



or 



dr 2 d0 



+ q rAj2, 1 = 0. 



Therefore 



or 



mr 



mr 



2 ^ 
dt 



dt 



+ q 

^ = 



= constant of the motion, 



constant/r • 



At r = 0 the left hand side is zero requiring that the constant 
be zero* Thus 



d0 

dt 



q A^/mx . 



(5) 
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Now substituting (5) into equations (1) and (3) gives 



d , dr, 
dt dt^ 



bA 



0 



m 0 Sr 



d , dz . 

5T <" dt> 



2 

a_ A 

m ^0 dz ' 



Upon rearranging the derivatives in these two equations one obtains 



d , dr . “Q * \ u x. ux r u. , . n 

dt^^ dt^ m ^0 Sr ^^dt^ 2 dz dt dz dt^ ^ 

dz 



2 ,2 .2 



0 _ fdz. d r dr r dz d 



dz , 



and 



d / dz. _ -q^ . ^^0 _ dz r d . dz . 
dt ^ dt^ m 0 Sz dt dz^^ dt^ ^ 



They can now be combined into the one equation 



m 



^ 2 ,2 , 2 hA^ 2 SAw 

i_^_^rS_A __^n+a_A —^=0 

^dt^ ^2 dz ^ m > 9z ' m ^0 9r 
dz 



( 6 ) 



dz 

The term (■^) is eliminated using the vector velocity relationship, 
equation (4), and equation (5), as follows: 



or 



dr ^ d0 ^ dz ^ 

dt dt ^ dt 



,d2,^ 2 ,di.^ 2 ,d0.^ 



and substituting gives 

d2,== 



(g) =[v" - (3 a^) 1/[ 1 Mf) 1 . 



dz 

Now replacing ('^) in equation (6) leads to the desired trajectory 



equation 



d^r 

dz^ 



1 



— [k^- 1 - — 

2 ^0 ' dl 



a A, 



0 



0 9z 



+ = (7) 



where K^= p^/q^. 
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Therefore one must choose the value of p/q to be used and know 

the relation for A = f(B^,z,r) developed in the previous section 

in order to integrate the trajectory equation. 

The trajectories were calculated by numerical methods on the 

available computer. The integration routine, originally written 
10 

by Gagliano , used a fourth order Runge-Kutta method. However, 
it was rewritten for the IBM 36O-67 computer using the Hamming's 
modified predictor -corrector method which is the IMPGC subroutine 
of the IBM System/360 Scientific Subroutine Package. The com- 
plete computer program is listed in Appendix IV. Figure 5 shows 
a set of trajectories calculated on the computer for ^Li^ + ^^Ne 
where the energy of the incident Li ion was 100 eV. As in the 
case of the uniform magnetic field there are two methods of de- 
tecting particles scattered at different angles: (1) by moving 

the detector and keeping the focusing magnet current at a specific 
value; (2) by keeping the distance between the scattering cell and 
the detector fixed while varying the current through the focusing 
magnet coil. Both methods were used and acceptable results were 
obtained from each method. 

5 . ANGULAR AND ENERGY RESOLUTION 
As is shown in the diagram (Figure 6), our experimental 
arrangement for measuring differential scattering cross sections, 
da (0), as a function of scattering angle allows a spread of 
angles AO at 0 to be admitted by the detector. Therefore there 
is some angular spread associated with every measurement which 
causes an uncertainty in the value of da (0) obtained from the 
measurement. It is seen from the sketch (Figure 6) that the 
angular spread depends on the size of AZ^, whose value is determined 
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Figure 5 . Numerically Computed Ion Trajectories for Li on Ne at 100 eV. 
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igure 6. Angular Spread Allowed by the Detector Apertures. 



by the detector apertures and 0. Thus, AO can be calculated from 

the relation between AO and AZ .. This relation is obtained from 

o 

- f(E,0,B,RSZ',C) (1) 



where E = energy of the incident ion 
0 = angle of scatter 

B = Magnetic field = f(I);I = coil current 
R* = target displacement from system axis 
= target displacement from system center 
C = a constant for each reaction including mass 
ratios, inelastic energy losses, etc. 



The coil current is constant during a given measur ement 'making 
dl = 0, therefore 



dB 




dl 



0 . 



Also, for a given reaction, C is a constant which makes dC = 0* 
Thus differentiating (1) leads to 



/az\ az 



d0+ 



fdz 

c 

3R' 



dR» + 



EOZ’ 



^3Zq 

V^^Veor' 



dZ^ 



Solving for the angular spread we get 

^az 



°/eR'Z* \ °/ER*Z» 



Se, 



dE+ 



'hz 

c_ 

9R » 



dR«' + 



ER*Z‘ 



EOZ' 



o 

az* 

/EOR’ 



dZ 



The values of the various quantities in the expression for 
the angular resolution must be calculated for each and every cross 
section that is measured. However as shown by example in Appendix 
I the terms arising from the energy spread and the finite target 
size increase dO by less than 10% for ion beam energy spreads 
less than 0,5 eV. Therefore the angular resolution is 



dO 




dZ 

o 



(3) 
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when the terras containing dE, dR*, and dz* are ignored. 

The second quantity of interest is the spread in energies, 
dE^, that particles scattered at 0 raay have and still be detected. 
As can be seen from Figure 6, this quantity depends upon AO at 0 
and also upon the energy of the incident ion, E, and the reaction 
parameter C. Therefore 

E = E (E,0,C) 
and upon differentiating 




If one considers only elastic processes, where the incident ion 
of raass in interacts with a target of mass M, then, due to the 
conservation of energy and raoraentum^^ , 



E = E 
s 



2 2 2 2 2 2 2 2 
M - m + 2ra cos 0 + 2m cos 0 HM - m + m cos 0 



(M + m ) " 



Using this 



raE \ 



and 



9E \ -2m E Sin Q 

s 1 _ s 

90 



E vL2 „2 . 2 ..2 






- m + m cos 0 



giving 




( 5 ) 
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Using equations (3) and (5) one can evaluate the angular resolution 
for a particular cross section. As an example, for the reaction 
^Li^ on He at 300 eV, I = 13-0 amps, dE = 0,25 eV and 0 = 39.6^, 



one obtains 








d0 


= 0 . 38 ° 


and 








dE 

s 


= 3.77 eV 



which are representative values of the reactions that were 
observed . 

6 . SOLID ANGLE 

The major advantage of this method of measuring differential 
scattering cross sections is the greatly enhanced solid angle 
obtained by collecting all particles scattered into a conical 
shell dO at 0. The solid angle for this system is 

d0 = 2TTsin 0 d0. 



As can be seen from Figure 6 all the particles scattered at 0 

must cross the axis in AZ at Z in order to be detected. This 

o o 

AZ^, then, determines d0 at 0 and, consequently, the value of 
the solid angle at 0. In the discussion of the angular resolution 
the expression 



d0 




dZ 

o 



was obtained. This is the relation used to calculate the solid 
angle. Thus 



dQ (0) = 2rr sin 0 



! S0 

[hZ 



dZ 
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However only 290° of the available 3^0° of the axial angle were 
open because of the way the scattering cell is constructed (see 
Chapter III, Section 5 , for a full explanation). Correcting for 
this and expressing in degrees/centimeters, the solid angle 



expression is written 



dn(o) = 0.0873 



sin 0 



do 

dz 



dz steradians. 
o 



This expression must of course be evaluated for every trajectory 
to be measured since it depends not only on the angle of scatter, 
0, but also on AZ and the derivative (dO/dZ All of these 

quantities depend upon the particular ion trajectory considered. 

Some representative values of the solid angle as a function 
of 0 at various focusing magnet currents are shown in Figure 7 . 
These curves are for the reaction ^Li^ + ^*^Ne at 100 eV incident 
ion beam energy. Notice the dependence of the solid angle on the 
angle of scatter and the focusing magnet current (magnetic field). 
When the magnet current is held constant the intercept distance 
changes for each angle of scatter causing the solid angle to vary. 
These values at the solid angle are typical for all reactions 
observed . 
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Figure 7* Solid Angle vs Angle of Scatter for Elastic Scattering 
of ^Li^ on Ne at 100 eV Incident Ion Energy 



CHAPTER III 



EXPERIMENTAL APPARATUS 
1. VACUUM SYSTEM 

The vacuum shell shown in Figure 8 consists of. aluminum 
cans with neoprene 0-rings at the connecting flanges, feed 
throughs, and end plates ^ 

Two six-inch diffusion pumps are used to maintain the 
vacuum* A National Research Corporation Model VHS6 diffusion 
pump having a 2400 liter per second pumping speed is connected 
to the detector end of the system. Water cooled baffles and a 
liquid nitrogen cold trap are located between the detector chamber 
and the diffusion pumpc The forepump is a CENCO HYVAC 45 mechan- 
ical pump with a 280 liter per minute capacity. A stacked half- 
moon baffle is placed in the foreline to reduce backstr earning 
from the forepump. At the opposite end of the system, where the 
ion source and the mass analyser are located, a Dresser Model 
DPD-6 diffusion pump is utilized. The Dresser pump has an l800 
liter per second pumping speed and is connected to a Welch Duo- 
Seal forepump. Water cooled baffles and liquid nitrogen trap are 
installed between the diffusion pump and the mass analyser chamber 

Dow Corning 704 diffusion pump oil is used in both diffusion pumps 

-7 

This arrangement provides ultimate pressures of 2.0 x 10 Torr 

-7 

in the detector chamber, and 1.0 x 10 Torr in the ion source 

chamber. With gas in the scattering cell at a pressure of 

-3 

2c 0 X 10 Torr the pumps maintain the detector chamber at 

2,0 X 10 Torr and the mass analyser chamber at 4.0 x 10 Torr. 

These pressures are low enough to insure that reactions outside 
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Figure 8. Experimental Apparatus 



4o 




the scattering cell are negligible- Tests also showed that the 
pressure in the detector chamber is low enough to. eliminate any 
pressure effects on the electron multiplier which served as the 
detector for the scattered particles. 

2, IW SOURCES 

In order to generate the ion beams used in the differential 
scattering cross section measurements^ two different types of ion 
sources were used: Lithium ions were produced by a thermal emitter; 
helium and hydrogen ions by an electron impact source. The per- 
tinent facts for both types of sources are reported here while a 
more detailed account of our developmental work on the electron 
impact source is presented in Appendix VI. 

Lithium Ion Source 

The lithium ion source used was essentially the source de- 

12 

scribed by Heinz and Reaves and is shown in Figures 9 and 10. 

It is a thermal emitter consisting of a very porous tungsten 
matrix into which the minerals beta -eucr ypt ite (Li^O'AI^O^ •2SiO^) 
or spodumene ^ ^ kl ^0 have been fused. The lithium 

ions are produced when the tungsten matrix is heated to tem- 
peratures in excess of IlOO^C by a separate filament (see Figure 
10),, The emitter was enclosed in a heat shield to keep filament 
power requirements low. The shielded emitter was mounted on a 
8-1/2 inch diameter stainless steel flange using a stainless steel 
mounting. In front of the emitter was an extraction electrode 
and a three electrode einzel lens. To extract the ion beam, the 
emitter, shield, and filament were placed at a potential corre- 
sponding to the desired beam energy. The filament power was 
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Figure 9 . Schematic of the Lithium Ion Source. 
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Figure 10. Lithium Ion Source Emitter Detail, 
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adjusted for an emitter temperature of about 1170^C which was 
found to be an optimum operating point. The effect of the emitter 
temperature on the total ion emission is shown in Figure 11. The 
extraction grid was grounded producing a beam of the desired energy 
The einzel lens apertures were lok centimeters in diameter and the 
lenses were 1*4 centimeters apart. The center lens was 4.0 centi- 
meters from the extraction electrode. These dimensions were ob- 
tained from the einzel lens relations given in Appendix VI and 

-7 

produced ion beams up to 1.0 x 10 amps xn the scattering cell 
when the emitter was operated at optimum output. The ion beam ob- 
tained from the source was found to have an energy spread of no 
more than 0.5 eV in the energy range 15 to 800 eV. 

This type of ion source has several advantages over plasma 
type sources: simplicity of operation^ small gas load on the system 
stable and highly monoener get ic beams of high purity. The major 
drawback of course is a limited number of species which can be 
produced by this technique^ 

Electron Impact Ion Source 

This ion source is one of a large family of sources which 
produces ion beams by creating a plasma in the appropriate gas 
and extracting ions from the plasma by externally applied elec- 
trostatic fields o In the particular source shown in Figure 12 j 
the plasma is produced by electron bombardment of the low pres- 
sure gas. Electrons j emitted from the filament, accelerated 
radially outward through the concentric grid and decelerated as 
they approach the cylindrical walls, oscillate between the fila- 
ment and the wall until they are collected on the grid. The 
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Figure 11. Lithium Ion Source Emission Curve. 
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Figure 12. Electron Impact Ion Source. 






plasma forms between the grid and filament and the ions are ex- 
tracted axially from it. To produce an ion beam, source gas is 
introduced into the source can at a rate that will maintain a 

-3 

pressure of about 1 x 10 Torr within the source. The potentials 
are set at the desired ion beam energy and the filament power is 
increased until enough electrons (usually about 90 nia ) are pro- 
duced to maintain the discharge. The grid is usually 100 volts 
above the common can and filament potentials. An einzel lens 

placed directly in front of the extraction hole helps produce 

-8 

ion beams up to 3»0 x 10 amps in the scattering cell. Our 

measurements indicate a rather large energy spread of about 8 eV 

although other workers using the same source have reported energy 

* 

spread as low as .5 eV . 

The source is fairly easy to operate and maintain, although 
requiring more maintainance than the lithium source. The primary 
drawbacks are the gas load it imposes on the vacuum system and the 
fairly large energy spread. 

A more detailed account of our developmental work on this 
ion source is presented in Appendix VI. 

3. MASS ANALYSER 

The ion sources used tend to produce a variety of ions- from 
the same source. This was especially true when using a plasma 
type ion source. Although the source parameters could be varied 
to enhance the production of one particular ion, other ions could 

*Private communication from Dr. W. Aberth, SRI. 
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not be reduced to insignificant values. Therefor e^ it was decided 
that mass a.''a lysis of the ion beam prior to its entering the scat- 
tering cell was necessary. 

The mass analyser shown in Figure 13 was designed and cali- 

13 o 

brated by Strohsahl . It has a 45 angle of bend in order to 
allow for both high beam transmission 3 up to and a resolution 

Am/m of about 40 allowing mass analysis of singly charged particles 
up to mass number 40 and energies up to one kilovolt. Using a coil 
current of 8 amps, it had a maximum field strength of about 4000 
gauss confined between 40 ^ wedge sector shaped pole pieces with a 
3/4 inch gap width. The pole pieces were external to the vacuum 
system which allowed increased flexibility and enhanced the vacuum 
integrity of the entire system. 

One effect resulting from the use of the mass analyser was 
the reduction of background signal at the detector by a factor 
of 100 or more. This fact made it worth while to use the mass 
analyser even in cases where the ion source did produce only one 
type of ion. 

4« FOCUSING MAGNET 

The magnet cell consists of five concentric brass spools 
wound with a total of 7938 turns of number 15 wire connected in 
such a way (see Table I) that the total resistance was approximately 
10c75 ohms. A single layer of 1/4 inch outside diameter copper 
tubing was placed over the outermost layer of wire on each spool. 
Water was run through the tubing for cooling. Thus it was possible 
to use up to 14 amps through the coils. A Sorenson Nobatron model 
OCR 15 O-I 5 power supply was used to supply the required magnet coil 
current regulated to + 0.1%. 
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SCALE INCHES 



Figure 13. Schematic of Mass Spectrometer. 
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TABLE I 



FOCUSING MAGNET DATA 



ITEM 

1. Number of layers of 
winding 

2. Average number of 
turns per layer* 

3 . Total number of 
turns on spool 
(item 1 X item 2) 

4 . Measured resistance 
( ohms ) 



ICOIL 
1 1 


COIL 

2 


COIL 

1 3 


COIL 

! A 


COIL 

5 i 


[total 


1 22 
! 


l8 


1 

1 16 


i 

12 

1 


lA ; 


82 


I 

j!97 

1 

1 


95 


97 


1 

|99 

1 


98 1 

1 

1 


1 

i 

1 


i 

2116 

i I 


1710 


1552 


1 

!ii88 

1 


; 1372 

j 


7938 

i 


i' 

! 20.46 


\ 

18.36 


20.86 


1 

118.04 


23.92 


10.75 



! 



5 . Length of wire (ft) i, 

(item 4 X 1000/2.58)** j 6010 

i 



6. Inside diameter of 
spool (inches) 



7.750 



6790 

11.155 



7560 '6790 9300 

14.300:17.250 19.750 



36 , 450 
7.750 



7. Inside diameter of 
winding (inches) 

8. Outside diameter of 
spool (inches) 



8.000 11.400 14.550 17.500 

i ! 

! 

11. i 4 o| 14.285 17.235 19.735 



20.000 

22.250 



8.000 

22.250 



*Based on actual count of several random layers. 

** Based on nominal resistance for 15 wire of 3.25 ohms/1000 ft 
at 25°C. 

coil 5 

I am 1 

coil 2 coil 3 

— oma iijisM— 

coil 1 coil 4 

^—WMJUL 

10.75 ohm 

COIL ARRANGEMENT 
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since the magnetic vector potential can be expressed in terms 
of the axial component of the magnetic field on the axis only this 

quantity was carefully measured. The measurements were carried 

7 . . . 

out by Kelly using a Hall probe calibrated in a magnetic field 

varied from 0 to 5 kilogauss against a known nuclear magnetic 
resonance probe accurate to 1 gauss at 1000 gauss. Figure l4 
illustrates the method used to position the magnetic field probe. 




TA-410522-II 



Figure l4. Arrangement for Measuring the Magnetic Field. 

The probe was placed inside an aluminum pipe on the tip of a 
wooden rule accurate to 0.5 mm. The wooden rule provided an 
adequate measurement of the distance of the Hall probe from the 
center of the magnet. The aluminum pipe holding the probe was 
fitted through templates at either end of a larger aluminum pipe 
held rigid inside the magnet coil and positioned concentric with 
the magnet coil axis. The end plates could be rotated to different 
radial angles and they had seven holes, 1.27 cm apart, along a 
radius which allowed field measurements at various r and 0 values. 
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The axial component of the magnetic field on the axis was 
measured at 4 centimeter intervals at current settings of 2,4,6, 
8,10,12, and l4 amps. A plot of these measurements is shown in 
Figure 15* These values were fitted to a 12 — order polynomial 
in Z and used to calculate the magnetic vector potential as ex- 
plained in Section 3 of Chapter II. Figure l6 shows the value 
of the axial component of the magnetic field measured at various 
r and 0 values at the center of the coil. These and other mea- 
surements tabulated in reference 7 demonstrated that the magnetic 
field was symmetric about both the center of the coil and the 
coil axis to within 0.1%. This allowed use of the coil axis 
and coil centerline as references for alignment of the entire 
system to the magnetic field axis. 

5o SCATTERING CELL 

The scattering cell, shown in Figures 17 and l8, constructed 
of stainless steel, is designed to allow particles scattered at 
all azimuthal angles to leave the scattering cell and eventually 
be collected at the detector* Figure l8 illustrates how this was 
accomplished by opening the scattered beam exit slit about 290 ° 
around the scattering cell. 

Attached to the front of the scattering cell is a circular 
aperture of variable diameter which can be adjusted from outside 
the vacuum system. Through appropriate wiring the ion currents 
to the front section of the scattering cell can be monitored. In 
the front section of the scattering cell a 0.5 centimeter diameter 
by 1.0 centimeter long beam entrance tube was drilled. These 
dimensions were chosen to reduce gas flow out from the scattering 
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AXIAL COMPONENT OF MAGNETIC FIELD Gauss 




-80 -60 -40 -20 0 20 40 60 70 

Z — cm 



TB-4I0522-I2 

Figure 15. Axial Magnetic Field Strength vs Distance. 
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Figure l6. Axial Magnetic Field Strength Measured at Various 
r and 0 Values. 
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Figure l8. Schematic of Scattering Cell and Beam Collector. 



cell and allow maximum flexibility in the 1.0 centimeter diameter 
by 1.0 centimeter long scattering chamber itself. Into the bottom 
of this chamber a 1/8 inch outside diameter copper tube was placed, 
through which target gas entered. A precision gas metering valve 
placed at the other end of the inlet tube was used to allow pre- 
cise control of the gas flow into the scattering chamber. A 
VG-IA ion gauge was located so that its entrance tube was placed 
90^ from the gas inlet tube. The narrowest portion of the path 
to the ion gauge was 0.5 centimeters in diameter, occur ing where 
it entered the scattering chamber itself. From the rear of the 
scattering chamber the front section of the scattering cell 
sloped upward to the outer edge of the cell at an angle of 49^ 
from the horizontal. This surface acts as the front edge of the 
scattered beam exit slit. The back edge of the slit was formed 
by the face of the rear of the scattering cell which sloped upward 
at an angle of 3^^. As shown in Appendix III these two angles 
determined the gas target thickness. The inside of the rear of 
the scattering cell was opened up to a 2,54 centimeter diameter 
in order to allow free travel of the unscattered ion beam to the 
beam collector. This after section is also wired so that ion 
currents to it can be measured. The rear section of the scat- 
tering cell can be moved horizontally from outside the vacuum 
system. In this way the target thickness could be varied. The 
opening could be adjusted to any desired length from fully closed 
to up to 0.8 centimeters when fully opened. 
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The front and rear sections of the scattering cell were placed 
inside a stainless steel tube ( 0 .l 8 centimeter wall thickness) and 
separated from it by teflon insulation. The support tube ran the 
length of the entire scattering cell but was cut away at the exit 
slits so that only 19 % of the tube was left to obstruct the scat- 
tered beam. This area was accounted for in the calculation of 
the solid angle. The tube has four variable length legs on each 
end (a total of eight). The individual legs in each set are 
placed 90^ apart around the circumference of the tube end. This 
arrangement gives very stable support to the cell and yet allows 
easy adjustment for alignment. 

The beam collector consists of three grids and a collector. 
Normally, the grid farthest from the collector was grounded, the 
middle grid was given s small positive potential to keep slow 
positive ions from reaching the collector, and the grid nearest 
the collector was given a small negative potential to suppress 
secondary electron emission from the collector. The collector 
was kept at ground potential. The beam collector and its grids 
are placed inside a metal cup which both shields them from ex- 
ternal charged particles and also acts as a gas cap over the rear 
end of the scattering cell. The beam collector assembly is 
mounted on a rod which allows the assembly to be pulled away from 
the rear of the scattering cell and then swung out of the path of 
the incident ion beam. This is done to allow alignment of the 
entire apparatus along the magnet axis and to allow the incident 
ion beam to fall directly on the detector for calibration purposes. 
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6. SCATTERED ION DETECTOR 



The differential scattering cross sections to be measured are 
extremely small, and therefore it will be necessary to use some 
sort of amplification so that the scattered particle current can 
be measured on an ammeter. A Model 3^6 Magnetic Electron Multi- 
plier manufactured by Bendix Corporation, Cincinnati, Ohio, was 
chosen as the beam detector and amplifier. This multiplier, de- 

l4 

veloped by Goodrich and Wiley , could provide reproducible cur - 
rent gains up to 10 , sufficient for the projected measurements. 

The multiplier, shown schematically in Figure 19, consists 
of two plane parallel glass surfaces each 5.5 centimeters long 
and placed 0.5 centimeters apart. The glass strips are coated 
on their facing sides with a thin, high resistance (approximately 

g 

10 ohms) layer of tin oxide and antimony. Attached to the top 
of the field strip is an entrance grid of etched nickel mesh, 1.9 
centimeters wide and 1.2 centimeters high, with greater than 90 % 
transparency. Attached to the top of the dynode strip and 
directly behind the entrance grid is a tungsten cathode 2.54 
centimeters long. Permanent Alnico V magnets placed outside both 
the field strip and the dynode strip produced a uniform magnetic 
field of 400 to 600 gauss in the space between the strips. This 
magnetic field is oriented perpendicular to the long dimension of 
the strips and parallel to their faces. At the bottom of the 
strips and between them a stainless steel anode is located. The 
entire structure is quite rugged and easy to handle. Its outside 
dimensions are 10.2 centimeters high, 3*35 centimeters long, and 
2.06 centimeters deep. 
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Anode 



Figure 19. Operational Schematic of the Bendix Model 306 
Magnetic Electron Multiplier, 
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In the multiplier's normal mode of operation a potential 
difference is applied the length of each strip in such a way that 
the end of the strips nearest the anode is the more positive end, 
and so that the field strip is positive with respect to the dynode 
strip. This arrangement, along with the magnetic field of the 
permanent magnets, produces crossed electric and magnetic fields 
in the region between the two surfaces. In this field configu- 
ration the motion of any electron released from rest at any point 
on the dynode strip is a cycloidal motion in an E x B direction 
(down the long dimension of the strip). Because the equipotent ial 
lines are slanted with respect to the dynode strip surface the 
electron strikes the surface before the completion of the first 
cycle. For sufficiently high applied potentials the electron 
has enough energy to cause secondary emission with a secondary 
emission ratio greater than unity. The secondary electrons thus 
generated repeat the process of the original electron, until most 
all electrons reach the anode and are collected. Therefore any 
current incident on the tungsten cathode is read out on the am- 
meter connected to the anode multiplied by some measurable gain. 

As the scattered ion detector in the experimental apparatus, 
the multiplier has to be operated in a magnetic field of up to 
200 gauss. The effect of this magnetic field causes the secondary 
electrons of the multiplier to drift perpendicular to the long 
axis of the dynode strip. Thus some of these electrons drift off 
the side of the dynode strip and never reach the anode, thereby 
reducing the gain of the multiplier. The rate of drift increases 
with increase in the magnetic field, causing the multiplier gain 
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to decrease with increasing field. Shielding the multiplier from 
the external field was attempted using M» -metal, netic, and conetic 
foils. This proved unsatisfactory for two reasons. First, the 
large amount of shielding material required so disturbed the 
focusing field that scattered particle trajectories were not pre- 
dictable. Second, the gain of the multiplier could not be re- 
produced whenever the focusing field was changed and then re- 
turned to the original value. This was probably due to the hy- 
steresis in the shielding material. The most satisfactory system 
has been to operate the multiplier inside an aluminum shield to 
protect it from random background and to choose the various strip 
potentials which give the highest gain and lowest noise. For a 
particular set of potentials the multiplier gain was measured at 
different focusing magnetic field values. This was done by varying 
the magnet coil current and keeping the multiplier at a fixed 
distance from the magnet center and also by fixing the coil current 
and varying the distance of the multiplier from the magnet center. 

both methods gave essentially the same gain for equal values of the 

1 1*1 

axial magnetic field. A plot of these measurements and a 10 — 
order fit of the data points is shown in Figure 20. The fitted 
value was used in the actual calculation of cross sections. An 
actual gain measurement is made as follows: with the multiplier 
installed inside the apparatus, the focusing magnet current and the 
detector distance are established at their desired values. Then 
the field strip and the anode of the multiplier are grounded while 
the cathode is connected to an ammeter. The ion beam from the 
source is directed onto the cathode. This current (which has to 
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multiplier settings 




Axial Magnetic Field — Gauss 

Figure 20. Multiplier Gain vs Axial Magnetic Field Strength. 
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be kept below 10 amps to prevent saturation of the multiplier 
output) is then measured. Then, without changing the current 
incident on the multiplier cathode the multiplier is placed in 
its normal operating condition and the anode current is determined. 
The ratio of the anode current to the ion beam current measured 
at the cathode is the gain. This process is then repeated for 
each new set of magnet currents and detector distances. 

Measurements showed that the gain did not vary with gas 

"4 

pressure up to 2.0 x 10 Torr , a pressure much higher than the 
multiplier ever encountered when data was being taken. The gain 

.6 .7 4 

was measured using Li , Li , He , and H ions and no ion or 

mass dependence was detected. The gain was found to be inde- 
pendent of the anode current up to 10 ^ amps. Above these cur- 
rents the multiplier tends to saturate since the anode current 
is approaching the strip current in magnitude. To prevent slow 
ions from entering the detector, a grid was placed over the de- 
tector shield aperture. Normally, a small positive potential 
(about 20 volts) was sufficient to screen out this background. 

The detector signal falls about 50% as the grid voltage is in- 
creased from 0 volts to about 10 volts. As the voltage is further 
incr eased, the detector signal remains essentially constant until 
the scattered ion beam energy is reached, then the detector sig- 
nal again decreases sharply. The grid also permits a crude anal- 
ysis of the scattered ion energy. 

After about two months operation the multiplier strips be- 
came coated with vacuum pump oil. This caused large increases in 
the multiplier noise resulting in false current measurements. 
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Cleaning the strips with a pencil eraser to remove the deposits 
reduced the noise and did not effect the multiplier gain. 

Subject to the limitations outlined above the magnetic elec- 
tron multiplier proved to be a reliable detector for the scat- 
tering apparatus. Although initial measurements of the gain were 
tedious, once they were made, subsequent checking showed the gain 
to be stable. 

7 . SYSTEM ALIGNMENT 

In order to detect particles scattered at a particular angle 
the detector must be placed on the focusing magnetic field axis 
at the point where the scattered particles focus on the axis. 

The program for calculating the intercept distance assumes that 
scattering took place on the axis at the center of the focusing 
magnet field. Therefore the accuracy of the angle of scatter 
determination depends on how well the scattering center and the 
detector are aligned to the magnetic field. 

Since the geometric centerline of the focusing magnet coil 
was found to be the magnetic field axis to within 0.1%, this 
geometric centerline was used as the reference for aligning all 
other components to the magnetic field axis. To perform the 
alignment, crosshairs were placed on either side of the magnet 
spool to mark its axis. Then a small laboratory laser was 
positioned so that its beam lay on the magnet axis established by 
the two crosshairs. The crosshairs were removed and the line 
established by the laser beam was used as the reference for 
aligning all other components to the focusing magnet axis. This 
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method made possible the aligning of the scattering cell, the 
detector, and the focusing lenses to the magnetic field axis to 
an accuracy of 1.0 milimeter radical displacement off the axis. 
Analysis of computed trajectories showed this error to be 
essentially undetectable. 

The incident ion beam was aligned by adjusting the position 
of the source and the mass analyser with the ion beam on until the 
maximum possible beam reached the scattering cell beam collector 
using no deflector voltages. 
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CHAPTER IV 



EXPERIMENTAL RESULTS 



1, PRELIMINARY TESTS 



The differential scattering cross section shown in equation 
(1) of Chapter I can be rewritten in terms of the quantities 
actually measured as follows: 



To measure a cross section, then, a target gas pressure and a 
target length are established in the system. The magnet current 
and/or the detector distance are varied to cover the range of 
scattering angles to be observed. These two quantities also 
determine the value of the multiplier gain and, along with the 
angle of scatter, the solid angle. Then at each setting of the 
magnet current and the detector distance, I_ and I. are observed. 
These values are used to calculate da(0) using the appropriate 
computer program listed in Appendix IV. 

In order to eliminate any spurious signals a series of mea- 
surements was undertaken to test the effect of the variables I., 
Psc ^ upon the measured value of da. For a particular angle 

of scatter 




2 



( 1 ) 



where 



I^ = current at MEM detector anode in amps 
G = MEM detector gain 

I^ = incident ion beam current in amps 
P = pressure of target gas in Torr 

P = 3*536 X 10 par t icles/cm'^-Tor r 
t = target length in cm 
dQ = solid angle 



da = K ^ ^ , where K = constant 

i sc 
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Now for constant incident ion beam current, I., and constant 
target length, t. 



I 



D 




( 2 ) 



The experimental data plotted in Figure 21 verify this linear 
relation. Notice that the data are linear but does not extra- 
polate to zero. The fact that does not go to zero is probably 

due to multiplier background noise. However this error is not 

-3 

important above 1.0 x 10 Torr pressure where the measurements 
were usually made. Another linear relationship which must exist 



IS 



^D = (c^ Ih 



(3) 



for constant target pressure and target thickness. Figure 22 
shows a plot of one set of experimental data which demonstrates 
this characteristic. Finally, if the incident ion beam current 
and the target pressure are kept constant then the detector 
current must vary linearly with the target thickness. This re- 
lation, which is 




is shown experimentally in Figure 23. These three curves show 

that our detector signal has the proper dependence on the various 

parameters and behaves like a scattered particle current. For the 

differential scattering cross sections measured the ion current at 

-7 

the beam collector was about 10 amps, the current to the rear 
of the scattering cell was about 10 amps, the ion current on the 
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Target Pressure (Psc) - 10 ^ Torr 

Figure 21. Detector Current vs Target Pressure 
at 40.42°. 



for 



6 .+ 

Li on He 




6 + 

Figure 22. Detector Current vs Ion Current fox Li on He 
at 40 . 96 °. 
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Detector Current (Iq)— 10 Amps 




Figure 23. Detector Current vs Target Length for Li on He 
at 39.1°. 
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front face of the scattering cell was about 10 amps, and the 
detector signal of the MEM output varied from 10 amps to 10 ^ 

amps. With no target gas in the scattering cell and essentially 
the same current to the beam collector the detector current drops 
to about 10 amps. Consequently the signal due to scattered 
ions was at least 100 times greater than the multiplier dark 
current . 

2. He"^ + He SCATTERING DATA 

Another important test of our technique is the comparison 
of data measured by our method to data obtained by more con- 
ventional methods. We therefore measured the elastic differ- 
ential scattering cross sections for He"^ on He at 300 eV from 
40 ^ to 52 ° and at 400 eV from 42^ to 54^, and compared our data 
to that taken by Lorents and Aberth^, as shown in Figures 24 and 
25 . The gross shape of our data and Lorents' data is essentially 
the same, but we do not observe the fine structure in as much 
detail as they do. This is due to two effects, both tending to 
increase the angular resolution (A0 — 1^ for our data, compared 
to aO — 0 . 5 ^ for the Lorents data). One effect is the increase 
in angular spread that occurs because of the anomalously large 
energy spread of about 8.0 eV in the incident ion beam energy. 

The other effect is the increase in angular spread brought about 
by the spiraling of the incident ions about the magnetic axis. 
This spiraling will allow the incident ions to enter the target 
area at small angles with respect to the system axis, thereby 
increasing the angular spread the detector sees. Nevertheless 
these results do show that our technique can give reliable cross 
section measurements on an absolute basis. 
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3. Li”^ + He SCATTERING DATA 

Finally, a measurement of the elastic differential scattering 

6 + 

cross section of Li on He was carried out* Since, for this 
system, the helium target atoms have a smaller mass than the 
incident lithium ions there are two possible energies at which 
the lithium ion may scatter at any allowable laboratory scat- 
tering angle* There also is a maximum angle of scatter in the 
'■{ 

laboratory system determined by the relation 

At this maximum scattering angle the laboratory elastic cross 
section becomes infinite^”^* This particular property made this 
system suitable for a check of the accuracy of our angular de- 
termination and system alignment. Since our experimental appa- 
ratus does not detect particles scattered at a particular angle, 
but rather at a particular vector momentum (recall that the tra- 
jectory is determined by the p x B forces), this peak in the 
laboratory cross section can be seen by setting the detector at 
a convenient distance from the scattering center and slowly de- 
creasing the magnet current. This allows us to measure the cross 
sections continuous^ from the high energy side to the low energy 
side of 0 without moving the detector out of the scattered 
particle beam. Thus we are able to determine closely the magnet 
current for the largest cross section observed and compare that 
magnet current value with current value for 0 calculated from 
the theoretical trajectories. The results of these measurements 
are shown in Figures 26 and 27 which are plots of the laboratory 
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Figure 27 . ®iab ^Li^ on He, = 300 eV. 
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differential scattering cross section vs the angle of scatter at 

200 eV and 3 OO eV respectively. Notice that in both cases the 

upper branch, corresponding to higher scattering energies, suddenly 

increases sharply as we approach 0 (4l.7°) and then decreases 

along the lower branch as the angle of scatter decreases. The 

fact that the cross section peak appears at the correct angle 

illustrates that our angular determination is correct and our 

system is properly aligned. The spread in the data is about a 

factor of 2 and is believed to be largely due to error in the 

multiplication factor G of the magnetic electron multiplier and 

will be discussed in more detail below. The same data is shown 

in the center of mass system in Figures 28 and 29* There appears 

to be some structure in both curves most easily seen in the smooth 

data in Figure 30, However, from this data, it is not possible 

to conclude whether this structure is due to instrumental effects 

or represents an oscillation in the cross section produced by 

quantum mechanical interference effects. For the homonuclear 

case (e,g, He^ + He) oscillations produced by the interference of 

l6 

gerade and ungerade states has been thoroughly studied . In 

+ 17 

heter onuclear systems (e,g. Li + He) Lichten has proposed that 
at very small internuclear separations a similar interference 
effect may occur since the active electron might find itself in 
a field resembling that of a homonuclear molecule. This, of 
course, would cause structure in the differential scattering cross 
section. On the other hand, there are reasons to believe that 
instrument effects cause the oscillations. In the 80 eV data, the 
large dip in the cross section at 95^ may be due to interaction 
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Extrapolated from Lorents and Aberth 
data taken at 0 < 90®. 
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Figure 29. da(0) vs 0 for Li on He, E = 120.0 eV. 
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Figure 30. Elastic Differential Scattering Cross Section for Li on He. 



of the scattered particle beam with the scattering cell exit slit. 

This problem can occur for particles scattered at less than 36 ^ 

in the lab« The dips at 125^ in the 80 eV data and at 130^ in 

the 120 eV data correspond to particles scattered at about 4l.7^ 

in the laboratory system. Consequently it is possible that the 

dips are due to inaccuracy in the calculated solid angle very 

near the laboratory cutoff angle. Notice however that regardless 

of the structure our data has slightly less slope and is above the 

line extrapolated from data taken at smaller angles by Lorents and 
2 

Aberth . Our data indicates that the center of mass cross section 
tends to level out as the center of mass scattering .angle increases. 
Thus there is essential agreement between our data at the larger 
angles and that of Lorents and Aberth at the lower angles. 

4 . EXPERIMENTAL ERRORS 

There are several sources of error in the measurements. One 
source of error which cannot be entirely eliminated arises from 
the spiraling of the incident ions about the focusing magnet axis. 

As already pointed out in Section 2 above this introduces an error 
in the angular measurement, and is most important in measurements 
where the cross section has large variations as a function of the 
scattering angle. However, the effect is minimized by proper 
alignment of the incident beam to the focusing magnet axis. The 
gain of the MEM, which is the scattered ion detector, is a rather 
sensitive function of the magnetic field strength of the focusing 
magnet. Thus any small error in setting the magnet current or the 
detector distance would result in actual gains slightly different 
from those used in the cross section calculations. The error in 
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setting the detector distance and the magnet current would also 
cause differences between the actually observed angle of scatter 
and solid angle and their expected values from which the cross 
sections are calculated. The detector distance can be set to 
+ 1 mm and the magnet current can be regulated to + 0.05 amp. 

For the ^Li^ on He data this results in a possible error in the 
scattered angle of less than 0 . 05 ^ and in the solid angle of less 
than 5%. Another source of error is in the measurement of the 
target gas pressure in the scattering cell. Although an accurate 
determination of this error was not possible it is apparent from 
the scattering cell geometry that the actual pressure is higher 
than the measured pressure used in the cross section calculation. 
It is felt that the primary source of scatter in our data is due 
to the multiplier which had to be operated at potentials different 
from the recommended values in order to obtain high gain in the 
magnetic field. An analysis of the cross sections for ^Li^ on He 
give a precision which is on the order of + 20% of the average 
value of a set of measurements. 
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CHAPTER V 



CALCULATION OF THE INTERATOMIC 



POTENTIAL OF THE Li - He SYSTEM 



+ 



1. INTRODUCTION 



The estimation of the interaction energy between two atoms 

as a function of their separation distance is of primary importance 

in the understanding of numerous chemical and physical problems. 

The interaction energy can be found in several ways; however the 

most direct method is to calculate it from elastic scattering 

cross sections. In this chapter the potential function for Li"^ 

on He was calculated from the cross sections measured by our 

2 

technique and some data taken by Aberth and Lorents . The pro- 
cedure used to calculate the potential function is based on 

3 l8 

Firsov^s method (outlined by Lane and Everhart ) using Gauss- 

Mehler quaduatures to perform the necessary integration. The 

potential function was then fitted to a Born-Mayer form and com- 

19 

pared to the potential function calculated by Zehr and Berry 
from their experimental data. 



relates the potential function, V(r), to the center of mass angle 
of scatter, Q, and the impact parameter, b, by the expression 



2 . THEORY 



Classical scattering theory, valid for this calculation 



11 



0(E,b) 




( 1 ) 



r 



o 



where E is the center of mass energy, 

r is the separation distance, and 
r^ is the distance of closest approach. 
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The necessary and sufficient conditions for scattering particles 

into the element of solid angle dQ = 2nsin0 dO is that the 

cm 

particles be incident in the annular ring of area 2TTbdbo There- 
fore the cross section cr(0) is related to b and Q by 



2ncr(0) sinO dO = - 2Tibdb 



or 



TT 



= 2\ cr 
0 



(Q) sinQ dO. 



( 2 ) 



Since the functional form of cr(0) isn’t known the integration 

cannot be carried out. However b can easily be found from a plot 

of cr(0) sinO vs 0. The area under the curve from rr to 0 is then 
2 

b /2. Thus the relation between Q and b can be obtained directly 
from the scattering data • However equation (1) cannot be used to 
obtain V(r) directly, so the following method will be used to 
determine the potential function from the 0 = f(b) result. Let 

2 V(r)-i 2 

q = [1 - E 

2 2 

V(r) = [1 - q A^] E (3) 



and by knowing q and r, V(r) is determined. Making this sub- 
stitution in (1) gives 

r -h 

0 = TT -J 2b [ (^)/dql [q^-b^l dq 
b 

00 

f 

= TT 2b (d Jin r/dq)(q^-b^) dq . 
b 
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Now using the fact that 



a 



2 , 2 , 



n b(dj0n q /dq)(q -b ) dq 

b 



leads to 



0(b) = bj [din(q^/r^)/dq^l[ q^- b^] dq^ 



where the integration variable has been changed to q . Now 

2 2^ 

multiply through by db/(b -q ) and integrate from q to infinity. 



Thus 



, 00 ,00 



= \ bdb \ 

q q b 



[dfn(qj^/r )/dq^] dq^ 

7^2 2 ,%, 2 ,2.h 

(b -q ) ^(q, -b ) ^ 



(4) 



The integral on the right can be transformed using Dirchlet’s 
formula to give 

(" Q(b)db _ 1 (” ^ C" d(b^) 

) .W 2 h 2 \ 2q 2^'^^l\ 2 2 A 2 2 h ‘ 

(b -q ) J 2 1 r ^q 2 (b -q ) (q^^-b 



The right hand integral can be placed in a standard form and 
integrated to give rr . Thus equation (4) becomes 

( 0(b)db n ( ^(%) ^,,2,2,, n „ 2,2 



Thus 



r(q) = q exp 



.00 

1 0(b)db 

TT 



(b^-q^)-^ 



(5) 
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and by evaluating this equation and using the results in equation 

(3) one can obtain V(r). The evaluation of the integral in 

21 

equation (5) is done using Gauss -Mehler quadratures , where 



s’ 



n/2 

f(y)(l-y^)"^ dy = ^ L f(cos^l^rr) 



and n = even integer . For our case let y = q/b, then 

/-I 



Q(b)db 

2 2 ^ 

(b^-q2)2 



t: - 



Q(q/y)dy 

y(i-y^)^ 



Thus 



f(y) = 



2is/Jl 

y 



3. PROCEDURE 

As an example of the procedure consider the calculation of 

6 + 

V(r) from the data for Li on He at 80,0 eV . The cross section 
is plotted in Figure 30* A linear fit of thisdata was used in the 
potential calculation. From this data a plot of a(0)sin0 is con- 
structed. This plot is shown in Figure 31 • The curve has to be 

extended to l80^. The area under the curve from any Q to l80^ 

2 

is b /2. The result of this integration is plotted as Q(b) vs b 
and is shown in Figure 32. From this graph r(q) is evaluated. 

A value of q is chosen and 0(q/y) is found for each term of the 
Gauss-Mehler quadrature formula. Then the individual f(y) terms 
are found and the integral is evaluated. From this one obtains 
r(q) and finally V(r) using equation (3)^ The Gauss-Mehler quad- 
ratures give the same answer to four significant figures for 

n = 8,10, or 12. However we used n = 12 for all calculations. 
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Figure 31. 



<J(0)Sin 0 vs 0 for 
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4 . RESULTS 



Using the procedure outlined above the potential function 

was calculated from our data at E = 80.0 eV and E = 120.0 eV 

cm cm 

and from Aberth and Lorents data at E = IO 9 . 1 eV and E = 72.7 eV 

cm cm 

The results are shown in Figure 33 where the potential is plotted 
against the internuclear separation. Notice that our data fits 
quite well with that of Aberth and Lorents. 

The plot of 109 ^ ^ found to approximate a 

straight line, thus the potential function could be represented 
by the Born -Mayer form as 

V(r) = 259e r in 10 ^ cm. 

7 + 

Zehr and Berry who also measured the scattering of Li on He 
found their data to fit a Born-Mayer potential which is expressed 

V(r) = 370e"^‘^^ 

This is plotted on the same graph with our potential (Figure 33). 

-8 

The agreement is quite good down to about r = 0.3 x 10 cm, 
below which their curve rises more sharply. This discrepancy at 
small values of the internuclear distance is not surprising in 
view of the fact that the measurement by Zehr and Berry was carried 
out by a difference method which is more accurate at the smaller 
angles of scattering corresponding to larger internuclear distances 
The agreement with the data of Lorents and Aberth appears to be 
excellent . 



89 



170 

160 

I 50 

140 

I 30 

I 20 

I I 0 

100 

90 

80 

70 

60 

50 

40 

30 

20 

10 

0 



Zehr and Berry 




\ 



a - 120.0 c V 
© - 109. I e V 
A - 80.0 t V 
+ - 72 .7 # V 




.2 



33. 



Intexatofnic Potential for Li^ -He as a Function of 
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CHAPTER VI 



CONCLUSIONS 

The most difficult and most time consuming portion of the 
work reported here was designing the apparatus and demonstrating 
the feasibility of using the focusing properties of an axially 
symmetric magnetic field to measure differential scattering cross 
sect ions • This was accomplished by showing that data taken by 
this method is in essential agreement with that measured by other 
means. Therefore this method with its much larger solid angle 
will allow future measurements of heretofore too-small-to-be- 
observed cross sections. 

The second accomplishment of this research was the mea- 
surement of the differential scattering cross section of ^Li^ 
on He at large angles, thereby extending the data about 80^ P^st 
previous measurements. These measurements also demonstrated the 
case in data collection that this technique (with its built in 
momentum analysis) offers for systems with heavy incident ions 
and lighter targets. Finally the interatomic potential function 
at distances less than iR for the Li^ on He system was calculated 
from our data. The agreement of this result with that derived 
from other scattering measurements is additional evidence for the 
soundness of our technique. 
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APPENDIX I 



THE EFFECT OF TARGET SIZE AND ION BEAM 



ENERGY SPREAD ON THE SOLID ANGLE 



The expression for the angular spread (see Section 5 



Chapter II' “ - 






dE 



dZ’ 



d0 = 



0R»Z' 



E0Z^ 



E0R » 




ER’Z' 



where E is the ion beam energy y 

dE is the energy spread of the beam, 
dR'is the effective target radius, 

dz*is the effective target length. 

I No analytical method exists for calculating the four partial 

derivatives. They are computed by solving a number of trajectories 

where E, 05 R*, and Z' are varied appropriately and the corresponding 

change in the intercept distance, Z^, is noted. The value of dE 

offers no problem since it is the energy spread of the incident 

ion beam. It depends only on the source being used and can be 

6 + 

measured independently (for the Li source dE — 0.25 eV). The 
quantities dR ’ and dZ’ cannot be arbitrarily chosen because they 
are functions of the trajectories and the detector geometry. The 
maximum size of the interaction region is determined by the dia- 
meter of the cell’s beam entrance aperture and the cell’s exit 
slit through which the scattered ions depart as illustrated in 
Figure l8 (Section 5 of Chapter 3). However, the actual region 
from which an ion can scatter and enter the detector is smaller 
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than this and depends upon the trajectory of the ion. This can 
be seen by recalling that the effective detector opening is de- 
termined by which depends on 0^ (see Appendix II), which in 

turn depends on the particular trajectory. Thus a particle can 

enter the detector only if it crosses the axis within Z + AZ /2 

0—0 

and at the proper angle, 0^, regardless of where it originated 
in the target gas. Therefore to evaluate dR » or dZ ’ one would 
have to evaluate trajectories for many different values of dR’ 
or dZ^ until the trajectory with the maximum dR ^ or dZ* is found. 
In general the derivatives (^Z^/^R*) and (5Z^/BZ*) are linear for 
a given reaction. Thus the change in dO due to varying values of 
dR’ and dZ’ is linear and symmetric about its value for dR* = 0 
and dZ’ = 0, The most extreme values of dR’ and dZ* an ion can 
have and still be detected will give the largest positive and 
negative corrections to dO. At values larger than these, dO is 
zero and, consequently, dQ is zero. Thus it is apparent that the 
finite target size corrections are such that the average dO is 
essentially equal to the value of dQ obtained for scattering from 
the target center. It is suggested that a detailed study of the 
problem be made on an analog computer by tracing the trajectories 
backwards from the detector to the scattering cell and defining 
dR’ and dZ*^ by this technique. The average solid angle is then 
obtained by integration over the total target volume contributing 
at each detector and magnet setting. 

However, one may obtain an estimate of the maximum effect of 
finite target size on the solid angle by calculating the change 
in AZ^ due to changes in the interaction volume, as follows: For 
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an ion to be detected, it must recross the axis within a length 

Z + AZ /2 where Z is the distance between the detector and the 
0—0 o 

scattering center and AZ is determined by 0 and the detector 

o o 

geometry (see Appendix II). Now let the detector be placed at 2^ 

so that it will detect particles scattered from the center of the 

cell at an angle 0. The angle the particles make with the axis 

as they recross it is 0^. Ions scattered at 0 at some point in 

the target other than the center will recross the axis at Z’ with 

o 

an angle 0^. Since AZ decreases as 0 increases those particles 
^ o o o ^ 

with 0*. > 0 will be detected only if 
o o 

AZ AZ 



Z 

o 



o 



2 



< Z’ < Z 
o o 



4 * 



O 



2 



and there will be no increase in the effective solid angle from 

them. However, those particles with 0’ < 0 will be detected when 
’ ^ o o 



Z 

o 




Z» ^ 
o 



Z 

o 



+ 



AZ^ 

o 



2 



and the effective solid angle will be increased. Therefore we 
need only to consider those trajectories where 0’ < 0 . An 
analysis of the trajectories of ions scattered from different 
places in the target show that those scattered on the ion source 
side of the center and those scattered off the axis will recross 
it at distances less than those scattered from the center and that 
they will have larger 0^ values. Consequently they cause no in- 
crease in AZ since they must cross the axis within AZ of Z to 
o o o 

be detected. Those scattered on the detector side of the center 
or those which originate off the axis and recross it before leaving 
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the scattering cell will cross at distances greater than those 



scattered from the center and will have a 

they will be detected if they have a ^ 

the actual AZ correction is 
o 



smaller 0 value. Thus 
AZ^ ^ 

Z + . Therefore 

O Z 



AZ 



c 

o 



AZ + AZ • 
o o 

2 



Therefore, for a particular energy, angle of scatter, and magnetic 
field one finds (30^/Sz^) for varying dR ’ and dZ’ values from 
trajectory results, then gets 



0» 

o 



Q + 
o 




) AZ 
^ o 



which gives a AZ^. This is applied to give the proper value of 

c c 6 + 

AZ . Table II shows the results of the AZ calculation for Li 
o o 

on at 100 eVo Notice that varying dR ’ gives the larger value 

for (BO /Bz ) and therefore the smallest value of 0». This was 
^ 00 ' o 

used to calculate AZ* since this was the largest correction. The 

o 

largest error is about 8% and as the angle increases the error 
decreases . 

The energy correction is straightforward. The derivative 
(Bz^/BE) is evaluated from trajectory data. The energy spread 
dE is determined separately for the particular ion source and 
the correction (Bz^/BE)dE is calculated. This correction is 
about 6%. 
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APPENDIX II 



CALCULATION OF THE INTERCEPT SPREAD 

The fact that the detector apertures are of a finite size 

allows all ions whose intercept values are 2 *+- AZ /2 to enter 

the detector and be collected. This allowable intercept spread, 

AZ^, depends on the detector aperture geometry and on the angle 

0^ between the ion path as it enters the detector apertures and 

the center line. This is readily seen in Figure A2-1. In 

general, then, as 0 increases, AZ decreases. 

^ ’ o ^ o 

The specific relationships between AZ and 0 can be found 
^ ^ o o 

as follows: 

Consider an incident beam of ions whose outer edge is 
determined by the outer aperture and whose inner edge is deter- 
mined by the multiplier grid opening. AZ^ for this case can be 
calculated from the geometry shown below. 




Now 

but 

and 



AZ 

o 

£ 

s+f 

2b 



( +a ) cotan 0 
^ ^ o 

s cotan 0 

o 



tan 0 

o 
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Figure A2-1, Detector Aperture Geometry. 
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or s = 2b tan 0 -f 

o 

giving H = 2b-f cotan 

which gives AZ = (f"*-a) cotan 0 -2b. (1) 

The maximum 0^ a particle can have is that angle defined by a line 
from the outer aperture edge to the farthest edge of the multiplier 
grid. For this trajectory AZ^ = 0. Thus 

AZ = 0 = (f+a) cotan 0 -2b 
o ' ' o 

cotan 0 = 2b/(f+a) 

o 

0 = cotan“^ ( ) 

putting the values shown in Figure A2-1 gives 0 , = 33-2^. 

^ ^ o(max I) 

As we decrease 0^ ^equation (1) holds until the angle of the 
incoming ion is the same as the angle formed by the edge of the 
outer aperture and edge of the shield aperture. This angle is 
found from the relation 

tan 0 = 

o b 

^ . -1 , f-h , 

0 = tan (—r—) 

o ' b ' 

and putting the values 

0 , . T\ ® / Tx\ = 26 . 6 ^. 

o(min I) o(max II) 

This is the smallest 0 for which equation (1) determines AZ . 

O V / Q 

Now consider an incident beam of ions with 0 = 26.6^ whose 

o 

outer edge is determined by the shield aperture and whose inner 
edge is determined by the multiplier grid opening. This is shown 
in the following figure. 
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Again 



and 



g iving 



A2 

0 

1 

h+s 

b 

s 



A2 

o 



- + a cotan 0 

o 



s cotan 0 

o 



tan Q 



o 



b tan 0 -h 
o 

(h+a) cotan 0 
^ ^ o 



b. 



( 2 ) 



As 0 decreases the limit on the inner edge of the incident ion 
o 

beam will finally be determined by the line from the edge of the 
shield aperture to the nearest edge of the multiplier grid. At 
this time 



tan 0 

o 



a-h 

b 



or ® / • TT\ “ ® / TTT\ “ 17-3 . 

o(min II) o(max III) 

This is the lower limit on the use of equation (2). 

The final case occurs when the incident beam of ions is 
limited only by the shield aperture. The following sketch 
illustrates these limits. 
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Here AZ = 2h cotan Q . (3) 

o o ^ ^ 

Thus, the three possible ways a beam of incident ions can be 
limited by our detector geometry along with the applicable ex- 
pression for AZ are: 

^ o 

I. For 33.2° s 0 ^ 26.6° 

o 

AZ = (f+a) cotan 0 -2b; 
o ^ ' o’ 

II. For 26.6° Sr Q s 17-3° 

o 

AZ = (h+a) cotan 0 -b; 
o ^ ^ o’ 

III. For 17.3° ^ 0 > 0° 

o 

AZ = 2h cotan 0 
o o 



Figure A2-2 shows a plot of AZ^ vs 
angular range* 



0 over 
o 



the appropriate 
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ViO- zv 




Figure A2-2. Intercept Spread vs Intercept Angle. 
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APPENDIX III 



CALCULATION OF TARGET THICKNESS 



The thickness, t, of the target depends on the scattering 
angle, 0, the slit openings, which can be controlled from outside 
the vacuum system, and the angles which the sides of the scattering 
cell make at the exit slit. 




a = 2.54 cm 
b = 2 . 00 cm 
h = 1 . 50 cm 



From the above sketch of the scattered particle beam exit, 
the four points (numbered 1,2,3>4,) that restrict the scattered 
beam can be seen. For particles scattered at angles less than 
36 ^, points 2 and 4 limit the scattering length. For particles 
scattered at angles equal to or greater than 49^j points 1 and 3 
limit the scattering length. 

Now consider t for particles scattered at angles, 0, less 
than 36 ^. Then t must be the distance along the cell axis between 
two parallel lines each making an angle 0 with the axis, where one 
of the parallel lines passes through point 2; the other through 
point 4> as shown on the following sketch. 
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To find t analytically, construct a third line parallel to the 
first two and require that it pass through point 1. 

Then t = s - f and I can be found from the law of sines 

^ _ a 

- sin 0 

sin (36 -Q) 

or jC = a s in ( 36 ^- 0 ) /s in 0 

and t = s - a s in ( 36 ^- 0 ) /sin 0. (1) 

Now consider t for particles scattered at angles 0 equal to 
or greater than 36 ^ but less than or equal to 49^. This case is 
shown in the sketch below. 
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Here t is the distance along the scattering cell axis between two 
parallel lines each making an angle 0 with the axis; one of the 
lines passing through point 1, the other through point 4. Then 
it is obvious that for this case 

t = S. (2) 

Finally, consider t for particles scattered at angles greater 
than 49^* For this case t must be the distance along the cell 
axis between two parallel lines each making an angle 0 with the 
axis; one of the lines passing through point 1, the other line 
passing through point 3, as shown below. 




It can be seen from the sketch that t = S - 
Again using the law of sines one has 

i ^ ^ 

Sin(Q"49^) SinQ 

giving f = b Sin( Q-49°)/SinO 

t = S -b Sin( Q-49^)/SinO. (3) 

In summary: Q < 36 ^: t = S -a Sin ( 36^-0)/SinQ 

36 ° ^ 0 ^ 49°: t = S 

49° ^ © : t = S -b Sin(0-49'^)/SinQ 
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A plot of vs 0 is shown in Figure A3-1. 



From this relation 



the scattering length, t, can be obtained. Figure A3-2 is a 
of t vs 0 for a scattering cell opening of O .3 centimeters, 
that it depends on two things: ( 1 ) the length of the scatter 

cell exit slit and, ( 2 ) the angle of scatter being observed. 



plot 

Note 

ing 
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Figure A3-1. Target Thickness Correction vs Scattering Angle (Lab). 
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Figure A3-2. Target Thickness vs Scattering Angle for an 0.4 cm Cell Opening. 



APPENDIX IV 



LISTING OF COMPUTER PROGRAMS 

The computer programs used to calculate the solid angles and 
the cross sections are listed in this appendix. All the programs 
were written in Fortran VI and were run on the IBM 36O-67 computer. 
The programs are: 

1 . SOLANG: This is used to integrate the trajectory 

equation and to calculate the solid angle. This 
program was written by N,R.A, Smyth. 

2 . CSVSZ: This program calculates the differential 

scattering cross section from the raw data where 
the observed angle is varied by moving the detector 
while keeping the magnet current constant, 

3. CSVSI : This program calculates the differential 

scattering cross section from the raw data where 
the observed angle of scatter is varied by changing 
the magnet current while keeping the detector fixed. 

4. PSIQ^A: This program calculates the differential 

scattering cross section from the raw data for the 
case where the target has less mass than the pro- 
jectile. 

All three cross section programs give both the laboratory system 
and center of mass system values for the cross section and 
scattering angle. 
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APPENDIX V 



ION BEAM OPTICS 
22 

Many good references on the design and use of electro- 
static and magnetic lenses to focus ion beam already exist and 
are well known by workers in the field. Therefore, only the 
space charge effect and the focusing properties of an einzel lens 
will be discussed. 

Space Charge 

The coulomb interaction among the charged particles in an 

ion beam causes the beam to spread as it moves away from the 

source. This spreading due to the space charge sets an upper 

limit on the amount of beam that can be transported between two 

apertures. The addition of converging lenses along the path of 

the ion beam will help reduce the space charge spreading and help 

deliver more beam to the target area. 

The maximum current that can be transported a distance x 

starting with a radius y and ending with a radius y can be 

c a 

calculated from the following assumptions: 

(1) The radial current density is constant across any 
beam cross section. 

(2) The current density is zero outside the outer ray. 

(3) The total current is constant. 

(4) The axial velocity of the beam is constant. 

The field E at the surface of the beam can be calculated from 
Gauss’ Theorem. Thus 

5pdv = fe E da 

^ y 
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wher e 



p = space charge density 



da = surface element 
dv = volume element. 



Now integrating at gives 

2 

p TT y Ax = C E 2TTy Ax 
c y 

since all the flux goes out through a cylindrical surface at 
radius y. Thus 

2 , 

E = Py /2ye 
y C 

2 

The current 1 through a cross section ny^ can be found from the 
linear velocity of the ions, u^. Therefore 



l/ny = Pu 

C X 



or 



P = lA^u y 
x^ C 



Now 



d^ 2 

qE = m 1 = q(py /2ye ) 

y dt^ ^ 



and substituting for p and solving gives 



d y ^ _3_ 



,.2 2mu 'rrey 
dt X 



(^) = b^/2y 



wher e 



b = ql/mu^nc = constant. 



Now 






dt^/ 



Ut = 






integrating 



dy 

dt 



= b^" log y + c. 



(1) 



133 



The initial radius was and at that time all ions were traveling 
parallel to the x axis, thus 



dy 



^ = 0 giving c = log y^. 



Substituting into equation (1) gives 



= b(log y/y^) ^ which is the differential equation 
for ions in the outer ray. To solve this equation let 
s = (log y/y^) 



or 



y = y^e 



Then 



2 

dy ,dy,,ds. , ^ s ds 

— — = (—^)( — ) = bs = V 2s e — 

dt Ms^Mt'' dt 



giving 



s^ ds 



Now 



2y. 



dt 



2 

{ s ds 
“ J ® dt 



or 






2y _ ' c 



s^ ds 

e -^ = D(s)-D(s^) 



dt 



23 

where D(s) is the Dawson function, which has been tabulated . But 

t-t = (x-x )/u , which 
c ' c' x’ 



gxves 



x-x 



(n —) b = D(S) 

^ 2y u ^ ^ ' 

C X 
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2u 



Thus 



(x-x )/y = -^ D(s). 



c' c 



2qV ^ 

Now using u = and also substituting for b gives 



m 



h 



4(ne) (i) 



(x-x )/y = oTTX 



Therefore the maximum ion beam current that can travel a distance 
starting at a radius and ending with a radius y, whose energy 
is qV, for ions having a charge q and mass m is 



3- Q/2 ^ 

I = l6ne (q/2m)2v^^ (D(s)) . 

c 

For our experimental apparatus 

x-x = 75*0 cm 
c 

y = 0 . 25 cm 
c 

y = 0.25 cm 



( 2 ) 



s = log y^/y = 1 and D(s) = 1.2463, giving 



I = 5.38 X 10 



-9 V 



3/2 



m 



where I is in amps, V is in volts and m is in amu. 

A plot of I vs V is shown in Figure A5-1 for several different 



ions . 

Einzel Lens 

In order to reduce the space charge spreading electrostatic 
lenses were placed at the source and also about half way between 
the mass analyser and the scattering cell. The three diaphram 
einzel lens shown in Figure A5-2 was chosen for both lenses. 
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Equipotential L»nes 
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Figure A5-2. Einzel Lens Equipotential Lines, 



In this lens the two outer diaphrams are set at ground potential 

and the potential on the inner diaphram is varied to give the 

maximum current at the scattering cell. Since the voltages on 

both sides of the lens are the same the refractive indices are 

equal on both sides and the three diaphrams act as a single lens 

with one focal length. In fact the focal length was found to be 

24 

independent of the applied potential. Chancon investigated this 
lens by ray tracing and found the following empirical formula for 
the focal length 



£= 6.5 

s 



S , -1/S. 

R <B> 

1 * I ) 



where s = lens separation 

R = center aperture radius 
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The distance between principal point P and mid-plane M is given by 

MP = 0.06 

o 

where = outer aperture radius 

T = thickness 6f the diaphrams • 

The focal length relation holds as long as T < S. A plot of 
the focal length relation is shown in Figure A5™2. For the 
lenses used in our apparatus 

S = 1.4 cm 

R = 0.7 cm 

which gives f = 6.3 cm 
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Figure A5-3- Einzel Lens Focal Length Relation. 



APPENDIX VI 



ELECTRON IMPACT ION SOURCE 

At the beginning of this experiment it was realized that 
several ion sources compatible with the experimental apparatus 
would be required. Therefore some extra effort was put into 
building an electron impact ion source and determining its pperating 
characteristics. Although the majority of the data presented in 
this thesis were obtained using a different ion source it is 
planned to use plasma type ion sources in future measurements and 
a summary of our developmental work is therefore included in this 
appendix . 

In this ion source, shown schematically in Figure a 6-1, the 
ions are produced by electron bombardment of the gas filling the 
source can. A plasma is formed in the'"-region between the grid 
and the filament and the ions are extracted axially from it by. 
the application of external electrostatic fields. 

Physical Description 

The ion source consisted of a 5.715 centimeter diameter in- 
verted glass press (see Figure a 6-1) made of number 3320 glass, 
with four, in line, equally spaced, l/l6 inch diameter pure tung- 
sten rods. The two outer rods support a 15 niil tungsten wire 
helical grid. The two inner rods support a 15 mil tungsten wire 
filament. Both grid and filament are 4.8l centimeters long and 
4.8l centimeters above the flat side of the glass press. In the 
original configuration the glass press became coated with a layer 
of tungsten after a very short period of operation. In order to 



Deflector Shims 




Figure a 6-1. Schematic of Electron Impact Source. 



I4l 



prevent the filament and grid from shorting, the glass press was 
inverted, letting the tungsten support rods sit in the well as 
shown in the figure^ With this arrangement no shorting of the 
source elements was experienced. 

Source gases were passed through a liquid nitrogen trap to 
remove water vapor before they entered the source. A precision 
gas metering valve placed in the inlet pipe was used to regulate 
the gas flow into the source. The source pressure was measured 
with a Bayar d-Alper t type ion gauge attached to the source by 
means of a 3/4 inch inside diameter ^^quick-connect” coupling 
attached to the source flange. This arrangement gave positive 
control and accurate measurement of the ion source pressure. 

The source can was made of number 304 stainless steel and 
insulated from its base plate with a 1.27 centimeter by 1.27 
centimeter lavite ring. The beam extraction hole was varied in 
size from 0.5 centimeter to 0.1 centimeter. An einzel lens con- 
sisting of three parallel stainless steel plates was attached 
directly in front of the ion source. The holes in the lenses were 
1.4 centimeters in diameter and the lens elements were 1.27 cen- 
timeters apart. These lens dimensions, calculated using the lens 
formulas given in Appendix V, were found to give the highest beam 
currents inside the scattering cell. Two pairs of steering elec- 
trodes are mounted external to the outermost focusing lens to 
assist in beam alignment. 

The ion source is wired as shown in Figure a 6-2. The fila- 
ment power is supplied through an isolation transformer which 
allowed the filament and the grid to float on top of the extraction 
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Figfure a 6-2. Electron Impact Ion Source Wiring Diagram^ 



electrode (can) potential. The highest potential is that of the 
grid and is the approximate beam energy. Energy analysis has 
shown the energy spread to be 4«0 eV about the grid potential, 
however other workers have found the energy spread to be + 0.5 
A filament current of approximately 11 amps was found to produce 
adequate electron emission of about 100 ma. This corresponds to 
a measured filament temperature of about 3800 ^K and a filament 
lifetime of over 100 hours. 

Operating Characteristics 

The gr id-cur r ent/gr id-voltage relationship, shown in Figure 
a 6-3 for 9 as, are characteristic of these of all gases used. 
The source operating area corresponds to the grid voltages above 
the grid current plateau. Although there is some measurable ion 
beam with a zero grid voltage, there appears to be a minimum grid 
voltage below which there is no actual discharge. These minima 
were found to be, for the different gases, 

GAS MIN. GRID VOLTAGE 

70 Volts 

He 90 Volts 

36 Volts 

-3 

at a source pressure gauge reading of 3«0 x 10 Torr » Any in- 
crease above these voltages caused a small, 1% at most, decrease 
in total beam intensity. 

The dependence of ion beam intensity on source gas pressure 
is shown in Figures a6-4 and a6-5 for hydrogen and nitrogen re- 
spectively. As can be seen from the figures, the maximum ion 
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Log Ion Beam Current in Amps Grid Current in Milliomps 





Figure a6-3. Electron Impact Source Characteristics. 
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Log (Ion Beom lnten§jty ) — Amps 




Log (Source Pressure Gauge Reading) Torr 

10 

Figure a 6-4. Ton Output vs Source Pressure for 
Source Gas. 
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Figure a6-5. 



Ion Output vs Source Pressure for N 2 
Source Gas. 
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beam intensity occurs at source ion gauge readings between 

-3 -3 

2 X 10 Torr and 3 x 10 Torr « Below this pressure range the 

ion beam intensity is quite small since there was no longer a 

discharge in the source. With the ion source operating properly, 

-8 

beam intensities on the order of 3-0 x 10 amps were obtained in 
the scattering cell. 

Theory of Operation 

The discharge in the source is a Townsend discharge, which 
is non-self-sustaining. There are two primary sources of ionizing 
electrons; one is the tungsten cathode which produces electrons 
thr.ough thermionic emission, the other is neutral atoms or mole- 
cules from which electrons are liberated in ionizing collisions. 
However, before these ^’secondary” electrons can ionize in turn 
they must be accelerated by the electric field until they obtain 
a kinetic energy equal to or greater than the ionization potential 
of the neutral gas. The effect of these secondary electrons can 
be seen in Figure a 6-3« There the large increase in grid current 
and beam current as the grid voltage is increased from 6o volts 
to 90 volts is explained by the fact that above about 60 volts 
the secondary electrons have enough energy to ionize the neutral 
gas through collisions. 

The beam current is also pressure dependent as is seen in 
Figures a 6-4 and a 6-5* This peak in the beam current vs pressure 
curve can be explained in the following way. An electron with a 
small energy in an electric field E begins a mean free path, L. 

To be capable of ionizing at its next collision it must gain energy 
from the electric field after traveling a distance d such that 
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qdE ^ q or d ^ — 



where V is the ionization potential of the neutral gas. The 
probability that the electron will travel a distance d without 
a collision is 



p(d) 



-d/L 

e • 



This is the number of ions produced by an electron in a mean free 
path. Since 1 cm of actual path contains 1/L mean free **ionizing” 
paths, the probability of a free path longer than L per cm of path 
in the field direction is equal to the number of ionizing col- 
lisions per cm of path in the field direction. Therefore the 
number of ions produced by an electron per centimeter is 

O' = P(d)/L = i 




The mean free path depends on the pressure and can be written as 
i = Ap , where A is a constant and p is the pressure. Then 

^ JL 

Ep 

Ae ^ 

a = p , where B = AV. 



From this expression, the pressure at which cv is a maximum, 
making the beam current a maximum, can be found. 



^ -B/Ep - ^ B -B/Ep ^ 

da=Ae ^dp-pA~e ^dp 

E 



(1 

dp 



E 
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At maximum 



- 0 making (1 - B/E p^) = 0 



Thus 



p = E/B 



is the pressure at which maximum beam current (ion formation) 
occurs. Notice that this depends on the gas, its ionizing 
potential, and the electric field the electron is in. This 
demonstrates that there should be a maximum in the pressure vs 
beam current data. However, our approach has been too sim- 
plified to allow us to calculate p^ from tabulated values of B. 
For example ionization by photons emitted from the tungsten 
filament has been ignored. However, from the simple model it 
is seen that in order to obtain the most intense beams, the 
source should be operated at or a little above p and with a 



grid voltage at or above that required to produce ionizing 
secondary ions. 

Proton Production 

Since a high intensity proton beam is desirable for the 

study of some of the proposed reactions an attempt was made to 

enhance the proton output of the electron impact ion source by 

introducing water vapor into the discharge as was suggested by 

25 

Dawson and Tickner The set of reactions 



m 





+ e 



»2° H 
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if actually occur ing, would enhance the amount of atomic hydrogen 
and hence the proton output while reducing the output. It was 
found that the proton output was increased with a corresponding 
reduction in the and outputs as water vapor was added to 
the discharge. However the total beam intensity itself de- 
creased by a factor of 100. It was also found that significant 
changes in the ratio of to and yields occur ed only after 
the discharge was 20% to 40% water vapor. It was felt that this 
was an unacceptable load on the vacuum pumps and this approach 
was therefore not pursued further. 
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